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1 Introduction 



Recent string theory computations [TJ [2] of the drag force on a heavy quark, and a related 
calculation [3] of momentum diffusion for non-relativistic heavy quarks, have raised the 
tantalizing prospect that the gauge-string duality might help us understand the dynamics 
of charm and bottom quarks propagating through the quark-gluon plasma (QGP) produced 
at the relativistic heavy ion collider (RHIC). Earlier work in a somewhat similar spirit 
includes [4]. Independently, a proposal was made in [5] for extracting the jet-quenching 
parameter q for M = 4 gauge theory from a Wilson loop calculation amenable to solution 
through techniques of classical string theory[j] Subsequent work includes extensions to non- 
conformal theories [HI [101 [IT], [12], non-zero chemical potentials [131 CHI EES! HSl EH HE], and 
other deformations of Af = 4 gauge theory [19]; studies of directional emission [201 HI] 
and the relation to the magnetic string tension [22]; and calculations of drag on particles 
carrying higher representations of the gauge group [23]. The venue for the string theory 
calculations is the gauge-string duality [24J [251 126] (for reviews see [271 [281 [29]), m particular 
the computation of Wilson loops through dual classical configurations of fundamental strings 
in anti-de Sitter space, as first considered in [3U1 l3"Tj. 

The drag force calculations [TJ [2J [3] result in a relaxation time tn for charm quarks that 
is in the right ballpark for comparison with data on the nuclear modification factor Raa 
and the elliptic flow parameter V2 for heavy quarks, as reported most recently by the STAR 
and PHENIX collaborations in [321 [33] . According to the prescription for comparing M = 4 
gauge theory and QCD advocated in [34], the string theory prediction is tu = 2.1 ± 1 fm/c 
for charm quarks. This value is lower than typical values ~ 4.5 fm/c used in one of the 
more successful phenomenological theories [35] of heavy-quark transport. But one should 
bear in mind that many steps separate the drag force calculation from predictions of Raa 
and v%. Prominent among these steps is a description of heavy quark propagation through 
the QGP which includes fluctuations in the force on the quark. The goal of this paper is to 
calculate the two-point functions of these fluctuations using the gauge-string duality. 

The organization of the rest of this paper is as follows. Section [2] sets the stage by 
reviewing the diffusion of momentum of a heavy quark in a Langevin formalism. Section [3] 
sets up the main calculation by finding the quadratic action and equations of motion for 
fluctuations of the trailing string. Section [4] shows how to translate appropriate solutions 
of these equations of motion into the Green's functions of interest. Section [5] shows how to 

1 There exists some debate in the literature regarding the calculation of [5]. Arguments from both sides 
of this debate can be found in [5J [71 [5] ■ 
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treat the zero-temperature case analytically. Section [6] presents a numerical study of the 
two-point function of primary interest. Section [7J includes a discussion of how the string 
theory results might be applied to understanding heavy quarks propagating through a real- 
world quark-gluon plasma. This discussion suffers from the usual difficulties of relating two 
significantly different theories, namely M = 4 super- Yang-Mills and QCD. Some significant 
technical issues are postponed to three appendices. 



2 Momentum diffusion in a Langevin description 

Following [36] one may attempt to describe the propagation of a heavy quark through a 
thermal bath in terms of Langevin dynamics (see also, for example, [3?1 138]): 

^i = - r]D (p)p i + F l L + F^ 

(if (* 2 )) = KMmWi-t*) (1) 

(if (ti)if (t 2 )) = «r(p) (6^ - p&Mh - t 2 ) , 

where pi = Pi/p is the unit 3- vector in the direction of the momentum p. The string theory 
prediction for the drag force (ignoring issues of modified dispersion relations raised in [1]) 
amounts to 

VD = ly/^N^, (2) 

independent of p, where m is the mass of the heavy quark. It is often assumed (modulo 
a subtlety having to do with how one discretizes the Langevin equations in the process of 
passing to a Fokker-Planck description) that kl and tjd are constrained by the Einstein 
relation 



k l = 2TE rj D = n^g 2 YM N T 3 7 , (3) 

where 

7 = 7rb (4) 

is the standard Lorentz factor for the heavy quark (and I have again ignored the possibility 
of a modified dispersion relation). The relation ([3]) emerges from requiring that a heavy 
quark propagating according to the equations (Cp) should eventually equilibrate to a thermal 
distribution e~ E ^ T , so it is a consistency condition for the Langevin approach. No such 
relation is available for Kt, but for non-relativistic quarks isotropy demands n T = n L . Indeed, 
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the result 



2T 2 / 

K L = K T = — = n^g 2 YM NT z fori;<l (5) 



was obtained through direct calculation in [3] and through use of the Einstein relation in [Tj. 

The Langevin approach (CQ) is hardly the only one in use for describing quark dynamics in 
the QGP. Another prominent paradigm hinges on radiative energy loss [39l HOI EFT] . Central 
interest in this description attaches to the BDMPS jet-quenching parameter q. The jet- 
quenching parameter as I prefer to define it is 

?T = (*>i)/A. (6) 

Here (p 2 j_) is the average transverse momentum acquired by a parton after it has traveled 
a distance A, measured in the rest frame of the plasma ^| The path length A is supposed 
to be taken large enough so that short-range correlations are washed out, but short enough 
so that the quark is only slightly deflected from its original trajectory. A definition of q 
that is preferred in some works, for example [5J [7j, refers to a partially light-like Wilson 
loop. To distinguish between the q of these works and the definition I prefer, I will employ 
a subscripted T as in ([6]). Note that the definition of qx does not require a strict light-like 
limit: it can be evaluated at any v. Also, does not require us to commit to the BDMPS 
formalism. For example, consider a Langevin description of momentum broadening where 
a heavy quark travels initially in the x 1 direction and is acted on by stochastic transverse 
forces satisfying 

d ^ = F l (FiitjFfa)) = S^Krit, - 1 2 ) , (7) 

where K^t) is an integrable function and i,j run over values 2, 3. If one assumes p±{0) = 0, 
then at sufficiently late times t (meaning times larger than the characteristic time-scales of 
Kxit)) one has 

<P±(t) 2 ) = S ij (p i (t)p j {t))= [ dt, [ dt 2 2K T (t 1 -t 2 )^2K T t, (8) 

Jo Jo 

where by definition 



k t = / dtK T (t). (9) 



2 1 would use pt in place of p± except for the fact that is usually reserved for momentum perpendicular 
to the beampipe. If the parton is traveling in the x 1 direction, then p± is the projection of the momentum 
onto the x 2 -x 3 plane. 
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Comparing with ([6]), one extracts 

q T = 2 -^. (10) 

V 

For a heavy quark with a known initial momentum, one may also consider 

q L = ((Ap L ) 2 )/\ (11) 

where Api is the deviation of the longitudinal momentum of a heavy quark from some 
average value. As with qr, the definition (TlTI) doesn't commit us to a specific formalism; 
however, a Langevin analysis precisely analogous to the one leading to (1101) gives 



Ql = ^ (12) 



V 

where 

/oo 
dt{F 1 (t)F 1 (0)) (13) 
-oo 

and iq is the fluctuating part of the longitudinal force. 

The AdS/CFT calculations that are the focus of sections [BE lead to force correlators 
from which qi and q\ may be extracted, essentially as in (Q, (TTOT) . (fT2j) . and ([TBI) . The 
results are 



g T = 2n^ g 2 YM NT 3 ^/v q L = n^/ g 2 YM N T^ 5 / 2 /v . (14) 

The result ( 1141) is larger by a factor 1/(1 — w 2 ) 3//4 than expected from ([3]). So Langevin 
dynamics does not capture the whole story. The method for obtaining (I14p is to numerically 
compute a symmetrized Wightman two-point function of oscillations of the trailing string 
solution of [TJ [2] and identify the zero-frequency component of this two-point function with 
kl or kt as appearing in (Q and ([TBI . The expressions in (fT~4"|) are obtained as a fit to 
the low-frequency limit of the numerics, but they are probably exact statements about the 
Green's functions in questionjf] 

3 Fluctuations of the trailing string 

The drag force prediction from the gauge-string duality, 



3 I thank C. Herzog for first suggesting to me the comparison with the non-relativistic results ([5]). 
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arises from a classical string configuration [TJ [2] in the AdSs-Schwarzschild background 

1 '-hdt 2 + dx 2 + 



ds l 



h 



7TT' 



(16) 



where T is the Hawking temperature. The string configuration, for p pointing in the x l 
direction, is described by 



x = x (t, y) = vt + 



(17) 



The fluctuating part of the force on the quark, which is located at the y = end of the string 
(i.e. on the boundary of ArfS'5-Schwarzschild) should be computable in terms of linearized 
fluctuations around the solution (117)1 . More specifically, if 



x 



1 ,1 



x\{t,y) + 5x 1 (t,y) x 2 = 5x 2 (t,y) x 3 = 5x 3 (t,y) , 



then the Nambu-Goto action 

>Sng = 



J d 2 a a/- det g a/3 g aj3 = G lu/ d a x' x dpx v 



2ixa 

may be expanded to quadratic order in 5x % to obtain 

L 2 1 



S 



NG 



2na' zh 



dtdy 



yj 



dtdyV a d a 5x 1 



l - I dtdy igfd a 5x 1 dp6x 1 + ^gfdaSx'dpSx*] , 

J V i=2,3 / 



where 



->a{3 4 potfi 



L 2 /z 2 H v (z H /y 2 (l-y 4 ) 



2na' y% y 

' z H 1 - yys 



L 2 /z 2 H 1 
2na' y 2 s 



> 2 \ 



y 2 (l-y*y 1-y 4 

2 44 
y y -y% 

y 2 z H J 



\ 1-y 4 



and I have introduced 



Vs 



(19) 



(20) 



(21) 



(22) 



(23) 
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If I switch from a a = (t,y) to some other parametrization of the worldsheet, then V a , Gl , 
and Q^P would transform not as tensors but as tensor densities: that is, they include a factor 
of a/- det g af3 . 

Because d a V a = and is proportional to Q^r , all three 5x l (t,y) obey the same 
equation of motion, namely 

d a (&? <f>) = (24) 
where <fi = 5x l , i = 1, 2, 3. Plugging an expansion 

(f>(t,y) = J ^0o(w)*(w,t,j/) V(u,t,y) = e- iut iJ;(u,y) (25) 

into (I24p leads to a "radial" equation for i/j: 

[s(y)d 2 y + t(oj, y)d y + u(u, y)] ^(u, y) = , (26) 

where 

t(u, y) = 2(1 - y 4 ) [1 - y 8 - v 2 (l - y" + zy 3 ^^)] (27) 
u(cj, y) = -yujz H [(1 - y 4 )uz H + v 2 y A (4iy + uz H )] . 

The radial equation (liZ6"j) has regular singular points at the zeroes of s(y): y — 0, y — x, 
and y = xys where x is any fourth root of unity The most interesting of these is y = ys, 
corresponding to some intermediate point on the string. What is special about this point is 
that it is the location of a horizon of the induced metric on the worldsheet. An intuitive way 
to see this is that a point on the trailing string with y held fixed follows a timelike trajectory if 
y < ys and a spacelike one if y > ys^ So the region y < ys corresponds to the exterior of the 
"black hole on the worldsheet," and y > ys corresponds to the interior. No signal from the 
interior can propagate classically to the exterior. Instead, signals in the interior region must 
by causality travel down the string. But there should be some Hawking radiation from the 
worldsheet horizon upward toward y = 0, and it is natural to suppose that it relates to the 
momentum diffusion. I will not make this connection directly, but the two-point functions 
that I will compute are related to both classical absorption and spontaneous emission by the 
worldsheet horizon. 

Although fl26l) does not seem to be explicitly solvable, certain limits of it are tractable. 

4 A closely related observation [TJ [42l [43j [6] is that there are no-drag configurations of mesons represented 
as strings with both ends attached to branes in ^ldS'5-Schwarzschild provided the string between the quarks 
hangs no lower than y = ys- 
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For example, to next-to-leading order in small y, two independent solutions are 



iMw, y) = 1 + + 0(y 4 ) </vev(^, y) = - ^|y 5 + 0(y 6 ) , (28) 



and to leading order in small positive ys — y, two independent solutions are 

$-{u),y) = 1 + 0(y s - y) M",v) = (.Vs ~ vf*"'** [1 + 0(y s - y)] . (29) 

Of the two solutions in ( l29l) . tf} + (u,y) corresponds to an outgoing wave, because the phase 
increases as one goes to smaller values of y. The standard horizon boundary condition, 
then, is to disallow this solution. To justify this statement completely, one should show that 
suppressing this solution amounts to stipulating that ty(u,t,y) should depend only on the 
infalling coordinate at the horizon. This is indeed true, because, to leading order close to 
the horizon, the infalling and outgoing coordinates are 

u- = t u + = * - \og(y s - y) . (30) 

A justification of (I3"U|) is postponed to appendix [Al 

Consider now a set of wave-functions ^n(u,t,y) = e~ luJt ipn(u,y), defined so that 

iPr(v, y) = ipDef (w, y) + C r (u)^ V ev(^, y) = (w)^-(w, y) ■ (31) 



In (j3T|) . Vr»ef) ^vev, and ip± are regarded as exact solutions of the radial equation (126]) . 
specified by their asymptotics at y = or y = ys- In the absence of an analytical solution 
to fl27|) . the quantities Cr{uj) and C^(u) must be determined through some approximation 
scheme or through numerics. The wave-functions ^>n(u,t,y) have several useful properties: 

1. ^r(lo, t, y) is a solution of the wave equation ff2^1) for 0. 

2. * fl (w,t,0) = e- <wt . 

3. \I/r(u;, i, ?/)* = to, t, y). To see this, note that complex conjugation is equivalent 
to sending uj — > —uj in the radial equation ff26l) . and that the boundary condition at 
the horizon is preserved by complex conjugation. This property guarantees that (p(t, y) 
is real everywhere provided O (^)* = 4>o{~ UJ )- 
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4 Momentum correlators from the trailing string 

As in the previous section, let <fi be one of the 5x l . Let O be the operator dual to 0, and let 
gaf3 _ gap QT ga/3 ^ ag a pp r0 p r j a t e- retarded, time-ordered, and symmetrized Wightman 
two-point functions are defined as 

G R (t) = -i9(t)([O(t),O(0)\) G F (t) = -i(T{O(t)O(0)}) 

1 (32) 
G(t) = -(O(t)O(0) + O(0)O(t)) 

where in general 

W; tre-^ ' 1 ' 

and H is the Hamiltonian of the gauge theory. 

In analogy to the recipe of [H] , the proposal for extracting the retarded Green's function 
from the wave-functions tyji(uj,t,y) is 

G R (cu) = -V R (cu, t, y)gy%y R (u, t, y) = log V R (cu, t, y) , (34) 



y=0 



where the notation \ y=0 means to evaluate at y = after taking the derivatives. The difference 
between (I3"l"|) and (3.15) of pH] is the use of rather than tpj^ Using ip^ would amount 
to restricting the sum over (3 to (3 = y, and this does not make sense in light of worldsheet 
reparametrization invariance. If the number current associated with ^u(uj,t,y) is defined as 

J a =^:G afl *R%*R, (35) 

then 

lmG R (u) = -J y , (36) 

and because J a is conserved, the right hand side of ( 1361) can be evaluated at any y. The 
number flux is positive at the horizon for uj > as a consequence of the infalling boundary 
conditions, so one finds from (|36|) that 

lmG R (oo)<0 forcu>0. (37) 

This is the correct sign for describing dissipative dynamics. If ipR had been used instead of 



5 Actually there seems to be one further difference: an overall sign. Possibly I have misunderstood the 
notation used [44]. Anyway, I claim that the sign in (f34|) is the right one. 
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tyji, the connection with the conserved number current would be broken. In [33] a restriction 
was made to consider only diagonal metrics, so the question of including dependence of the 
wave-function on coordinates other than the radial one never arose. (The context was slightly 
different: instead of Q al3 , the metric of interest was the line element of the bulk spacetime.) 
But (1341) is clearly in the spirit of [33] and its antecedent [25]: note for instance the similarity 
with (27) of [25]. 

The quadratic part of the on-shell action reduces to 

Son-shcii = \J M-^)Mu) ReG R (u) , (38) 

where to obtain the right hand side from the quadratic part of Sng I used an integration 
by parts, the equations of motion, the formula (1341) . and the property G R {— uj) = G R (u)*, 
which is obvious from property [3] at the end of section [3j 

Thus, just as in the examples treated in [33], the real part of G R (ui) may be extracted 
using a conventional AdS/CFT formulation in which the on-shell action is the generating 
functional for Green's functions. But it turns out that the imaginary part will be of greater 
interest, due to its connection with the dissipative dynamics. The relation between lmG R (uj) 
and G(uj), 

G {u) = - coth Tff— Im G R (u) , (39) 



is modified from the usual one in that Ty s appears in place of T. A justification of fl39|) . 
following [35] , is postponed to Appendix iBl 

Here is a somewhat informal argument for identifying 0(t) with minus the force on 
the quark in the direction that selects. To be precise, let <fi = 5x 2 : then the claim is 
0{t) = —p 2 (t). Consider the action 



S q [<Kt)] = J dtL^o) (40) 
for an external quark propagating through the thermal medium along a path 

x 1 (t)=vt x 2 (t) = 0oW- (41) 

I have excluded explicit dependence on O an d t from L because the medium is assumed to 
be translationally invariant and static. Note that translational invariance does not imply 
conservation of transverse momentum: S q is evaluated under the path integral and includes 
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couplings to the gauge theory. For small deviations 0o from a straight path, 

s q~ J dt ^°^~ = J dt ^ v ^ = ~ J dt 0oP<l> ■ ( 42 ) 

Of course, p^ = p2- Altogether, the generating functional for the two-point functions of 
interest is 

Z[(f> ] = (exp S^-iJ dtfop^J , (43) 

whence the desired conclusion O = —p\. More formally, p-iii) should be regarded as the 
operator which generates a local displacement at time t of the path of the Wilson loop for 
the external quark. A similar line of argument applies when = Sx 1 , and of course Sx 3 is 
equivalent to Sx 2 . 

Other Green's functions than G r {oj) could be computed by using other wave-functions 
than ^R(u,t,y) in fl34"|) . provided they satisfy the three properties at the end of section [3j 
For example, wave-functions ^a{^, t, y) = e~ tuJt ipA{^,y) with ipA proportional to ip+ should 
lead to the advanced Green's function. 



5 Zero-temperature Green's functions 

As a check on the validity of the prescription (I3"4"|) . and as a way of understand the u ^> T 
behavior in the T^O case, it is useful to detour to a consideration of the zero temperature 
limit, where the background metric is just 

ds 2 = ^-(-dt 2 + dx 2 + dz 2 ) (44) 

and the unperturbed solution is a string that hangs straight down: 

x 1 = vt x 2 = x 3 = . (45) 

I will use worldsheet coordinates a a = (t,z). An expansion entirely analogous to fl20]) can 
be performed, only now the integration measure is dtdz, and 

v * =0 <=d-^f=^(; _; /72 ). (4 6) 
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The wave-functions satisfying infalling boundary conditions at the degenerate horizon at 
z = oo were found analytically for v = in and can be adapted immediately to the case 
of interest: 

V R (u, t, z) = (! _ ilu)Z ) . (47) 

The retarded Green's functions are 

GJM = ic«M = (48) 

where I have dropped a divergent term proportional to uo 2 whose Fourier transform is a 
contact term supported at t — 0. The associated symmetrized Wightman functions are 

G t (uj) = -(sgnu)ImG*(u) = G L {u) = P~;\uj\ 3 . (49) 

lixa lixa 

In preparation for finding the real-time forms Gr(^) arid Gl(£), consider the general power- 
law Fourier integral: 

dt^= cZt^P = 2|a;| 2A - 1 r(l-2A) S in7rA. (50) 

■oo 1^1 J — oo I ^ I 

Expanding around A = 2 gives 

oo e iu;t 

dt — — = — \uj\ 3 + (analytic in to) . (51) 
t 6 

There is an analytic continuation implicit in the result (|5T|) : whereas the integral (150!) can 
be evaluated for < A < 1/2 with no regularization, the integral in (15T1) is highly divergent. 
Regularization schemes differ in how they prescribe the analytic terms on the right hand 
side of f l5TT) . corresponding to contact terms in real time, supported at t — 0. Employing 
f l5T|) and ignoring contact terms, one obtains 

Gr(t) = - 2 G L (t) = 3 -^^. (52) 

To me the most surprising feature of the final result (fl4|) in the case of non-zero tem- 
perature is the strong dependence of the magnitude of longitudinal fluctuations on 7. One 
can see already in (HHj) and (|52l a very strong dependence: oc 7 4 for fixed u; or t. In 
this zero-temperature setting, the factors of 7 should be capable of being understood just in 
terms of Lorentz invariance. As a step in this direction, consider the following toy model: 
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M = 4 super- Yang-Mills with gauge group U(l). The action for the heavy quark is 

S = J ds [-m\x\ + + |x|# 7 X 7 ] , (53) 

where dots denote d/ds and the parameter s need not be the proper time r. The 8 1 form 
a unit vector in R 6 , and they specify properties of the quark: its scalar "charges." For an 
external quark, one specifies the 8 1 as functions of s along with the trajectory The 
only case I will consider (which is also the case considered in [TJ [3], [2] ) is the one where the 
9 1 are constant. The classical equation of motion following from (1531) is 

m«„ = W M + e'X 1 ^ , (54) 

where I have now specialized to the case where s is the proper time r, and 

BX 1 

The velocity vector uF" = dx^/dr satisfies 

V uV = _1 Viw = diag{-l, 1,1,1}. (56) 

Let's now promote F^it^x) and X 7 (t,x) to operator-valued fields and consider the sym- 
metrized Wightman two-point function of W M (r). Because of conformal invariance and the 
obvious identity W^u 11 = 0, the answer can only be 

G™{t) = i(W M (r)W,(0) + W,(0)W m (t)> = C w ^±^ (57) 

for some constant Cyy- 

Whether W M is the entire "force" acting on the quark depends on whether the last term 
in Q54p is considered a force or a modification of the momentum. This is immaterial for 
external quarks because their mass is infinite and their trajectory can be prescribed to have 
constant . So for the purposes of comparing (|57|) to f[5"2"j) . this term doesn't matter. What 
does matter is that, in deriving (|52|) . the force was defined in terms of a t derivative, not a 
r derivative: 

F to = ii2, w (58) 

at 7 dr 

Using (jBTj) and fl58|) . and assuming as usual that the quark moves in the x l direction, one 
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obtains 



G L {t) = ^(Fi(t)Fi(O) + F 1 (0)F l (t)) = ±G™{th) = 7 4 ^ , (59) 



and similarly Gx{t) = ^Cyy/t 4 . Thus there is perfect agreement with (1521) if one formally 
makes the replacement Cw — * 3L 2 /7r 2 a'. 

The simple-minded analysis leading to fl59l) of course makes no prediction about the 
scaling of the low- frequency parts of Gl(uj) and Gt{u) with 7 when the temperature is 
non-zero. But what drives the relation Gl(u) = ^ 2 Gt{oj) is the corresponding relation 
Q2 = 7 2 ^^, and this holds whether T is zero or not. 

6 Numerical results 

Starting from fl22|) . fl3Tl) . (pi]l . and fl39l) . a short calculation shows that 

G*{u) = -Y(2C R (u) + iX) G T {uo) = Y coth (2 Im C^) + X) (60) 



where 
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i = /tV y = ^ — t— = a — z^9 2 y M N. (61) 

The longitudinal Green's functions are always 7 2 times the transverse ones. Only Gt{oj) will 
be evaluated explicitly here. To do so, it helps to note that the substitutions 



(62) 



/ 1 1 \ iw«£r/4 

v 2/ <*>*h 
y = — to = 

lead to a simplified form of the equation of motion fl26|) : 

[y (1 - y 4 )<9 2 - 2(1 + y 4 - ^y 3 )dy + w 2 y] x (w, y) = o . (63) 

An infalling solution ipn(ou,y), satisfying the asymptotics (13T1) . corresponds to a solution 
Xr(w, Y) satisfying 

w 2 Y 2 

XR ( Wy Y) = l + — — + c R (w)Y 3 + 0(Y 4 ) for small Y 

2 (64) 

XR ( W , Y) = c H (w) [1 + 0(1 - Y)} for y close to 1. 
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From (1601) . (|62|) . and (|64|) it is straightforward to show that 



and that 



where 



g *m - ( 2csW + 1™ 1 (65) 



Gt( U ) = 1 ^5 9T ( ro ) (66) 



Z H 


ys 


7TW\ 






jTW =(=coth^)(l + i^^M). (67) 

To calculate cr(«;), and thereby gx{w), I employed numerical integration, implemented with 
Mathematica's NDSolve, matched to high-order series expansions for Y close to and 1. 
Asymptotic forms at large and small w were extracted by fitting to the numerics: 

gr(w) = 1 + 1.26w; 2 + 0{w A ) for small w 

i, „/ 1 _„.\ . . (68) 

Graphs of gr{w) and its Fourier transform, 



gr{w) = — |u>| 3 I 1 + 0(w 8 ) ) for large w. 



° dw iwi. 



g T (£) = / e -^g T (w) , (69) 

J-oo 27T 

are shown in figure [TJ By visual inspection, the scale I* for the characteristic structures in 
g T {£) is t « 1.5. 

Because of the leading |w;| 3 behavior in gx{w) at large u>, the real space form <?t(^) is 
highly divergent at small £. This leading divergence is precisely the behavior observed in the 
zero-temperature Green's function in section [51 Contact terms supported at £ = restore 
integrability at the origin. A simple strategy for finding gr(^) is to subtract the |u;| 3 term 
from gr{w) and perform a numerical Fourier transform on the remainder, cut off at some 
maximum frequency W. The trouble with this strategy is that the subleading l/|w| behavior 
in gr{w) leads to errors of order log(W / |£|) at small £. Subtracting this subleading behavior 
as well leads to a remainder which is absolutely integrable at large w but not at small w. 
The simplest form I could find that agrees with the large w asymptotics in fl68l) to the order 
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9t(w) g T (i) 




Figure 1: The thick red lines come from numerical evaluations of gr{w) and the most accurate 
form of gr{£) I was able to obtain. The thin black lines are the analytic approximations 
g%? pTOX (w) and g^ pprox (£). The thin dash-dot line in the gr(£) plot is the residual, grity — 
gj? pTOX (£). The thick dashed line is an approximation to gr(^) obtained by subtracting off only 
the leading |u>| 3 behavior from gr{w) and numerically Fourier transforming the remainder 
out to W = 20. 

shown and leads to an everywhere absolutely integrable remainder is 

«r-W = lH"^5(i-^) (?o) 

9T PPrOX W = Y^"' (l + 4) "here I,, = (3/8)'/". 

The intuition of Debye screening suggests that gr{£) decays exponentially at large £, similarly 
to g^ pmx (£). 

As explained in [H], singularities of the retarded Green's function in the complex u 
plane correspond to quasi-normal modes of the corresponding field on the gravity side of 
the AdS/CFT duality. In the present context, the quasi-normal modes in question are 
fluctuations which are purely infalling at the worldsheet horizon and behave as y 3 for small 
y. Such fluctuations occur at values of uj where Cr(u) diverges. Equivalently, in terms 
of w = uzfj/ys, the quasi-normal modes arise at poles of cr(w). A preliminary numerical 
analysis yielded the following approximate values of the first two quasi-normal modes w n \ 

W! = 2.620 - 2.302i w 2 = 4.650 - 4.316i . (71) 

For each w n , —w* n is also a quasi-normal frequency. (Quasi-normal modes of the trailing 
string were also considered in [1] in the v — > limit, although the boundary conditions were 
different.) A standard expectation is that l\ = 1/Imtoi ~ 0.44 sets a scale for exponential 
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decay of the retarded Green's function. This is considerably smaller than note however 
that i\ relates only to asymptotic statements at large £, whereas is an approximate char- 
acterization of the main structures visible in grity- (Also note that because | Imw!| > 2, the 
singularity of gr{w) in the complex w plane closest to the real axis is not Wi but ±2z coming 
from the factor coth(7rw;/2), so a decay gr(^) ~ e -2 '^ is expected rather than e - '^ 1 .) 

7 Discussion 

One should be able to extract the transport coefficients kl, k>t, qr, and qi from the low- 
frequency limit of the symmetrized Wightman functions: 

v 

k t = ~qr = lim G T (w) k l = vq L = lim G L (u) . (72) 

The logic supporting f!72|) is a treatment of momentum diffusion as in equations fl7|)- (|T0|) . 
only with stochastic averages replaced by the symmetrized Wightman functions. Using fl66l) 
and \im w ^ Q g T (w) = 1, one finds the results quoted in ([TJJ: 

qr = ^ = 2n^NT^/v q L = ^ = n^NT^/v . (73) 

Because of the dramatic failure of the Einstein relation ([3]), it doesn't make sense to plug 
the trailing string predictions for rj^, k^, and Kt into a Langevin description of a finite 
mass quark: it won't equilibrate to a Maxwell-Boltzman distribution, due to the largeness 
of kl as compared to r]o at highly relativistic speeds. Perhaps what is needed is a stochastic 
treatment of the quasi-normal modes of the trailing string, the lowest of whose frequencies 
were found in (I7TT) . 

Because strongly coupled M = 4 gauge theory is different in significant respects from 
QCD, there is considerable uncertainty in how one should translate results like (T5]) and f!72l) 
into quantitative predictions for QCD. To characterize this uncertainty, consider the following 
two comparison schemes: 

• An "obvious" comparison scheme is to equate the temperature and the Yang-Mills 
coupling. I will assume T = 250 MeV and a s = 0.5, which are in a representative 
range for RHIC physics. To summarize: 

"obvious" scheme: T N=A = T QCD = 250 MeV g^ M N = 12na s = 6vr . (74) 
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• An "alternative" comparison scheme was proposed in [34] . The first part of this pre- 
scription is to equate energy density instead of temperature. Thus in any formula (for 
example (I73p ) in M = 4 into which T enters, one eliminates it in favor of a power 
of the energy density before comparing to QCD. This is approximately equivalent to 
setting T/v =4 = Tqcd/3 1 ' 4 . The second part of the prescription is to determine <?y M iV 
by matching string theory results for the force between a static quark and anti-quark 
to lattice results. This match is conspicuously imperfect, but at Tqcd = 250 MeV a 
range 3.5 < gyM^ < 8 was found by comparing static potentials at r pa 0.25 fm, with 
9ym-N = 5.5 suggested as a typical value. To summarize: 

"alternative" scheme: 3 1/4 T jV - =4 = T QCD = 250 MeV g YM N = 5.5 . (75) 

Combining (1731) with either (lT4"j) or (ITS"]) leads to the following estimates of qx for a charm 
quark propagating through a real QGP: 

"obvious:" q T w (2.2 GeV 2 /fm) 1 = 5.9 GeV 2 /fm for p c = 10 GeV/c 

vy 1 — u 2 

"alternative:" g T « (0.51 GeV 2 /fm) 1 = 1.4 GeV 2 /fm for p c = 10 GeV/c, 

1 — v 2 

(76) 

and for q L : 

"obvious:" q L w (1.1 GeV 2 /fm) _^ 2 ^ /4 = 150 GeV 2 /fm for p c = 10 GeV/c 

"alternative:" q L w (0.26 GeV 2 /fm) — ^ /4 = 36 GeV 2 /fm for p c = 10 GeV/c, 

(77) 

where I have set m c = 1.4 GeV. The choice p c = lOGeV/c corresponds roughly to = 
5 GeV for non-photonic electrons. The values of qx and scale with temperature as 
(Tqcd/250 MeV) 3 in both the "obvious" and "alternative" prescriptions. Based on the large 
values for qi in (i77|) . especially using the "obvious" comparison prescription, a conservative 
conclusion is simply that longitudinal fluctuations in the force are strong and merit a more 
thorough investigation. 

There appears to be some tension between the results of the present work and those 
of [SI tZ]. We insert a color source in M = 4 gauge theory in the same way, namely by 
letting a fundamental string end on the boundary. This is why our results have the same 



V 9ymN scaling. But the functional dependence as well as the magnitude of our results 
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differ significantly: in [5j [7J it was found for partons in the adjoint representation that 

vr 3 / 2 r(3/4) 



?lrw = r(5 ^ } ' y/9 2 YM NT* (71 



in the light-like limit, v — ► 1. Combining (1751) with either (I74p or (1751) leads to 

"obvious:" g LRW « 2.6 GeV 2 /fm 

(79) 

"alternative:" g LRW « 0.61 GeV 2 /fm . 

Again it is worth noting the sensitive dependence (Tqcd/250 MeV) 3 of ^lrw- It is also 
important to recall that ^lrw is defined differently from c[t, in terms of a partially light-like 
Wilson loop rather than momentum broadening. 

As a consistency check on the comparison of qr to a transport quantity used in modeling 
a real-world quark-gluon plasma, one should inquire whether the path lengths of real charm 
quarks are longer than the typical time scale t* of force-force correlators. An estimate of 
this time scale can be read off from the quantity discussed following (169]) as = £*zh/vs- 
Combining this with either (1741) or (175]) leads to 

"obvious:" U = (1.0fm/c)V#/10GeV 

, (80) 

"alternative:" U = (1.3 fm/c) tJE/10 GeV . 

Clearly these times are short compared to the lifetime of the QGP, which is roughly 6 fm/c in 
central collisions. But it is more relevant to compare to the mean path length of the charm 
quarks which escape from the QGP and are detected. This is presumably bounded below 
by the corresponding quantity for hard partons without a flavor tag, estimated in [47] as 
Amean = 2.5 fm. Perhaps an even more relevant comparison is the relaxation time to = 
for charm quarks as calculated via the trailing string picture. This quantity was estimated in 
[34] as to = 0.6 fm/c using the "obvious" prescription and 2.1 fm/c using the "alternative" 
prescription. The comparison with to is the most stringent: 

"obvious:" — = 1.7y/E/10 GeV 

't (81) 
"alternative:" — = 0.63^/10 GeV . 

tD 

The conclusion from (15T1) is that, depending somewhat on assumptions, time correlations in 
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the fluctuating forces on a charm quark probably matter, especially at higher energies. It 
would be useful to inquire what happens when the quark mass is made finite explicitly in 
the string construction by inserting a flavor brane as in [HI [T]. 

For b quarks, whose mass is roughly 4.8 GeV, time correlations are less important for the 
energies attainable at RHIC: t* scales as while tjj scales as m, so for fixed energy, the 

ratio t*/tf) is decreased by a factor of (mf,/m c ) 3 / 2 « 6.3. 

Correlation times may also help clarify the relation with [5j [7]: diverges as v — > 1, so 
if the path length A is held fixed as this limit is taken, it conflicts strongly with the u — ► 
limit in f!72|) . It is tempting to speculate that somehow cutting off the integral (Q at finite 
time would allow one to define a g-like quantity with a finite v — > 1 limit. 

A characteristic feature of radiative descriptions of energy loss following [39l |4T)| |4"T] is 
the scaling of AE not as path length A but as A 2 . The Brodsky-Hoyer inequality [4"9] . 



derived from uncertainty principle considerations, is saturated to within a factor of order 
unity in the BDMPS formalism. This is in contrast to a Langevin description including 
drag force, where AE oc A. As a rough comparison between a radiative description and 
the trailing string description, it is interesting to ask above what critical path-length A c the 
inequality (1821) is satisfied if one extracts the left-hand side from the drag force (ignoring 
longitudinal fluctuations in the force) and sets q = q^. The answer is 



For p c = 10 GeV, one finds A c = 1.0 fm using the "obvious" scheme and A c = 1.4 fm using the 
"alternative" scheme. In principle one can distinguish between predictions of radiative and 
Langevin descriptions of energy loss by making independent measurements of energy loss 
and momentum broadening. It has already been suggested [50] in the context of a Langevin 
description that two-point correlators in azimuthal angle for D-D pairs could provide some 
insight into the drag force. Further investigation of the potential uses of heavy quark two- 
point correlators to distinguish between competing models of energy loss would clearly be of 
interest . 

The trailing string makes the distinctive prediction that, for heavy quarks, both the 
average drag force and fluctuations are enhanced by powers of the Lorentz factor 1 / yl — v 2 
for v close to 1. This clearly implies that heavy quarks feel dissipative effects more strongly 



AE< l -q\\ 
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when they are more energetic. So the expectation from the trailing string picture is that 
Raa as a function of pp for non-photonic electrons has a more negative slope than predicted 
by treatments such as (5Tj in which q is held fixed (i.e. independent of pr)- It is therefore 
gratifying to observe in figure 3 of [33] a persistently negative slope of Raa over a wide 
range of momenta, 1.5GeV/c < p^ < 7GeV/c. The predictions of [51] include a weaker 
dependence of Raa on p^ over this ranged) The predictions of [SB] are also rather flat in this 
interval, despite a dependence of rjoip) on p that is only slowly decreasing in the range of 
interest. The predictions of [35], which, like |38j . are based on a Langevin approach, have 
a significantly larger momentum range in which Raa is strongly decreasing, and ?7dG°) as 
presented in the closely related work [37] is again slowly decreasing in the range of interest. 
Only three points, with p^ roughly between 4.8 and 6.5 GeV/c, lie below the predictions of 
[35] ; the same three, plus perhaps one other at pr ~ 4.2 GeV/c, lie below the predictions of 
[5T] . In summary, the comparison of Raa as reported in [33] with existing phenomenological 
models reinforces the view that dissipative effects increase with the momentum of the heavy 
quark; but without some more systematic study, it is hard to tell whether the trailing string 
picture does a better job of explaining the data than other approaches. 

Results from the STAR collaboration [32] agree fairly well with [51] if the contribution 
from bottom quarks is removed. But in light of the comparisons of pp — > eX PHENIX data 
[52] with fixed-order-plus- next-to-leading-log calculations [53], it seems ad hoc to neglect 
bottom contributions completely. The authors of [51] cogently warn of the difficulty of 
disentangling the contributions of charm and bottom in the absence of vertex reconstruction. 
In any case, the STAR data are at least consistent with a more negative slope for Raa as a 
function of Pt than either the c + b or c only curves from [51] . 



8 Version history 

In the initial version of this work it was claimed incorrectly that qr scales as \j\J\ — v 2 
as v — > 1. The correct scaling — v 2 was worked out in private communications with 

J. Casalderrey-Solana and D. Teaney, who independently obtained key aspects of the current 
work. Subsequently and independently, they and I worked out the case of longitudinal 
fluctuations. The crucial correction to my original treatment of transverse fluctuations is the 
replacement coth ^ — > coth in (1311 . as justified in appendix IBI following [45] . Inspired 

6 In referring to the predictions of one or another theoretical study, I am relying upon the portrayal of 
those predictions in [33] or [32] ■ This means, in particular, a choice q = 14GeV 2 /fm in the calculations of 

15a. 
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by their treatment, I realized that the Green's functions could be expressed as functions 
of w = uozh/us U P to overall powers of zjj and ys- This facilitated an analysis of the 
characteristic time-scale of force correlations. It also occurred to me that a simple zero- 
temperature treatment sufficed to understand the leading behavior at large w of g<r(w). 
Finally, in private communications, H. Liu, K. Rajagopal, and U. Wiedemann all urged me 
to distinguish even more carefully between q as computed through a Wilson loop formalism 
and qr = Taken together, these developments seemed an appropriate occasion for 

the substantial rewrite that occurred between arXiv versions 2 and 3. 



Acknowledgments 

I thank C. Herzog, H. Liu, J. Maldacena, G. Michalogiorgakis, C. Nayak, K. Rajagopal, 
U. Wiedemann, and especially J. Casalderrey-Solana and D. Teaney for useful discussions. 
I am grateful to B. Zajc for pointing out reference [50] to me, for correspondence regarding 
the electron's typical momentum fraction in D-meson decays, and for comments regarding 
the interesting dependence of Raa for non-photonic electrons as reported in [53]. I am 
indebted to J. Friess, G. Michalogiorgakis, and S. Pufu for collaboration on work that led to 
the results appearing in appendix El This work was supported in part by the Department 
of Energy under Grant No. DE-FG02-91ER40671, and by the Sloan Foundation. 



21 



A Infalling and outgoing coordinates 



A crucial property of the wave-functions n(u,t,y) appearing in (125|) is that they depend 
only on the infalling coordinate close to the horizon. The purpose of this appendix is to 
demonstrate that the infalling and outgoing coordinates are quoted correctly in (1301) . to 
leading order in small ys — y, where as usual ys = \/l — v 2 . Actually, I will do somewhat 
more by explicitly constructing coordinates in which the worldsheet metric is seen to be 
locally conformal to R 1 ' 1 (as any two-dimensional metric with — h signature must be). 

It is convenient to start by parametrizing the worldsheet using x l and y rather than t 
and y. The worldsheet metric is easily seen to be 



ds\ = g a (3da a da^ 



z 2 H y 2 



(dx 



1\2 



2W 



dx 1 dy + z 2 u dy A 



(84) 



where as usual h = 1 — y A . The manipulations to cast the metric ( 184|) in a conformally flat 
form are relatively simple because g a p depends on y but not x 1 : 



L 2 1 



y A /y 2 

1 - h/v< 



dy* 



1/ 



z 2 H y 2 



-Aidx 1 



;dy 2 



y 2 h — v 2 
n 2 {-dq 2 + drf) 



L 2 

o 2 J h - v 2 )dq 2 
v z H y A 



ZHyVV ;dy 2 



h/v< 



(85) 



where 



and I have defined new coordinates 



L 2 h~v 2 
y 2 v 2 z 2 H 



56) 



q = x 



dy 



zny 2 v 
h — v 2 



V 



±vz H Vl — v 2 J 



dy 



h-v 2 



57) 



Because h — v 2 = y% — y 4 , the integrals in (1H71) diverge at the horizon. So one must define q 
and T) piecewise in two regions: 



z H v I 1-iy/ys , ., l + y/ys\^ 
qi = x + — — log — h i log 



4iy s \ ^- + w/ys 



i-y/ys 



zhv 
Aiy s 



, l-iy/ys , ., i + y/ys 

lOg ; h % log 



i + iy/yt 



i - y/ys 



> 



for < y < y s ; 



(88) 
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X + 



4% 



log 



log 



i + iy/ys 
w/ys 



i log 



z log 
1 + 



1 + 2/A/S 

-1 + 

y/ys 



for y s <y <l, 



(89) 



V"° 1 + fy/vs -l+y/y S/ 
The second form of the metric in (1851) is analogous to the standard form of the Schwarzschild 
metric. In region I (outside the horizon) y is the spacelike variable and q is the timelike 
variable. In region II (inside the horizon) y is timelike and q is spacelike. 

Using x 1 = vt + £(?/), one may extract expressions for qi in terms of t and y rather than 
x 1 and i/; 771, on the other hand, is a function only of y. To construct infalling and outgoing 
coordinates in region I, one simply forms 



q± = qi± m 



(90) 



Expanding in small y$ — y, one finds 



vt + ... 



q+ = v 



t - 7T- log(ys - y) 
2y s 



+ ... 



(91) 



where the omitted terms comprise constant terms (i.e. independent of t and y) and terms 
involving positive powers of ys — y. Up to the overall factor of v, the expressions in (19T!) are 
precisely the ones given in (1301 ). 



B Relations between real-time thermal propagators 

The usual relations among the Green's functions defined in ( 1321) may be summarized as 
follows^] 



G f (lu) = Re G r (lu) + % coth Im G R {to) 
G{u) = -lmG F {u) = -coth^lmG R {u) 



(92) 



The first equality in (19*21) usually can be derived by inserting complete sets of states; the 
second is immediate from the definitions; and the third is an obvious corollary of the first. 
However, the usual derivation of the first equality relies crucially on having a fully equili- 
brated system, where (Q) = tr(e~ l3H Q)/ tre~^ H for any operator Q, and the heavy quark 
is far from equilibrated. It turns out that a very simple modification of (1921) suffices: one 



7 I am indebted to J. Casalderrey-Solana and D. Teaney for pointing out that the second equality in 
equation (25) of the original version of this paper was in error, and for private communications related to 
the contents of this appendix. The treatment in this appendix follows in its essential points that of |45| . 
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must make the replacement T — > Ty$ in (192]) . where ys = vl — v 2 . This applies equally for 
transverse and longitudinal fluctuations, and I will not need to distinguish between them in 
this appendix. 

A derivation of the modified form of fl92|) may be carried out using the Schwinger-Keldysh 
formalism, as adapted to AdS / CFT in [51] following earlier work [551 E3 E7] . For simplicity, 
let's set Zh — 1, so T — 1/ir. First define Kruskal coordinates U and V such that 

UV = y^± e ~ 2tan ~ ly - = -e u . (93) 
y + 1 U 

To make the most efficient use of these relations, it helps to note the following points: 

1. One may extend the second relation in (193]) over the entire Penrose diagram by as- 
signing t an imaginary part which is piecewise constant but changes by an odd-integer 
multiple of in/ 4 as U or V switches sign. A consistent way of doing so is exhibited in 
figure [2K 

2. The first relation in (1931) defines a monotonically increasing map from UV G (—1, e _7r ) 
to y G (0, oo), and UV = maps to y — 1. 

3. As a result of the previous two points, the relations (I9"3"|) define a one-to-one map 
(U, V) — > (t, y), where the domain of (U, V) is points in R 2 subject only to the restric- 
tion — 1 < UV < e _7r , and the range of (t,y) has four disjoint pieces corresponding to 
the four regions (RFLP) in figure [2k- The R piece, for example, comprises points in 
R 2 subject only to the restriction < y < 1. 

The trailing string solution (fT7j) may be recast in the form 

1) 

x 1 = - log V + v tan" 1 y . (94) 

(Recall that zh = 1.) In this form it is clear that the trailing string is entirely non-singular 
at the horizon between the R and F regions. The singularity at V = 0, corresponding to 
Ret — > — oo, is genuine. In the spirit of the Schwinger-Keldysh formalism, one prepares a 
state of the system at Ret = — oo and propagates it first forward along Imt = and then 
back along Imt = —tt/2. Correspondingly, one may use (194]) in the R and F regions and 

x 1 = -log(-V) + v tan' 1 y -ivrr/ 2 (95) 

in the L and P regions. The two solutions (1MI) and (195]) can be matched at V = by 
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1 □ 0<y<£ 

2 d e<y<y s -e 

3 ESI y s - £ < y < y s 

3 y s <y<y s +e 

4 EH y s + £<y<l-£ 

5 @ 1 -£<y < 1 

5 M l-£<y<l + £ 

6 D 1 + £<y 



(A) 



(B) 



Figure 2: (A) The Penrose diagram for AdSs-Schwarzschild with conventions on the phase 
of t specified. Increasing the real part of t corresponds to moving upward on the R boundary 
and downward on the L boundary. Increasing y means moving to the left in R and to 
the right in L. (B) The odd-numbered regions are narrow slices of the Penrose diagram 
corresponding to series expansions around y = 0, ys, and 1. The even-numbered regions fill 
in between these series expansions. Note that 3 and 3 are disjoint, but 5 is a subset of 5. The 
corresponding regions —1 to —6 are the reflections of 1 to 6 through the point U = V = 0. 



continuing V through the lower half-plane — that is, for Iml^ slightly negative. Such a 
continuation is consistent with the conventions laid out in figure [2] because in the L region 
one may rewrite fl95l) in terms of t and y in precisely the original form (fT7|) . 

Let's now consider perturbation wave-functions ty a for oscillations of the string in a 
slightly more general fashion than in section [3J The index a will run over values in a rather 
large set: 

a e {±1, ±2, ±3, ±3, ±4, ±5, ±5, ±6} . (96) 

These values correspond to different regions of the Penrose diagram, as indicated in figure [2b. 
In each of the odd-numbered regions, which are all narrow slices near some definite value 
of y, one may expand the general solution to the wave equation for \1/ to leading order as 
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follows: 
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h Bxy 3 ]e 
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(97) 

For brevity I have dropped subleading terms from the A-type solutions in ([97)1 even when 
they are bigger than the S-type solutions. There is no loss of generality in assuming e~ luJt 
dependence because the entire spacetime has a timelike Killing vector which is d/dt in each 
patch. 

The expressions in flDTj) involve imaginary powers of quantities which are real and positive. 
Such expressions are unambiguously defined: x lot = e ialogx . To specify ^ globally one must 
prescribe via some analytic continuation how the pairs (\l/3, ^3), (\l/-3, ^-3), and (^g, ^ _g) 
fit together. Also one must solve the wave equation for \l/ in the intermediate regions ±2, 
±4, either through some approximation or via numericsjf The result of such prescriptions 
and approximate solutions can be cast in terms of 2 x 2 matrices: 

v a = M Q/ 3 v/3 where v a = a I . (98) 



Here a and j3 are in the set ( 1961) . excluding the even numbered values. Note that M Q/ 3 for 
each a and f3 is a 2 x 2 matrix: a and j3 label the matrix, not its components. An obvious 
property of the M Q/3 is = M~^. 

The properties of the matrices M Q/ 3 that I will require are: 

M_ lj _ 3 = M 13 M_ g> _ 5 = M 35 (99) 

"^"33 I „ _ m ,,/o,,„ I "^"—3,-3 I „ —,.1/0., „ I (1^^) 





8 One could also find numerical solutions in region ±6 and perhaps match to some expansion near the 
future and past spacelike singularities. This would perhaps be of interest in an exploration of the physics 
of the black hole singularity, and it could also be of some relevance in computing higher point functions of 
operators pertaining to the heavy quark, but such matters are beyond the scope of the current work. 
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Mg i _g = e~/ 2 P JJ (102) 
M 13 D = ( 1 1 | = Mi 3 , (103) 

\c R H c R { u y-ix) 

where D is some diagonal matrix, Cr(u) is the same quantity as introduced in fl28l) . and 
X = 2v 2 uz H /3(l - v 2 ) as in (JHU). I will justify (|99D- fTT03|) after showing how to obtain the 
first equation in (l§2"j) from them. 

To calculate Schwinger-Keldysh propagators G Q p(u)), with a and j3 equal to ±1, one 
should specify A\(uS) and A-\{u) and calculate the "response" functions, Bi(u) and B-\{yj). 
To this end one requires the matrix 

M i.: = M i. M .'!. 3 M„3,„ 5 M ,. , M §,§ M-, , M 53 M 33 M 31 

~ A \~ (104) 
= Mi 3 (M - 1 ) 2 M31 = M 13 M^ 1 . 

Si ' \ o e ™/vs J 

The second equality in (11041) uses (}99l) . fllOOj) . ( 1101 h . and (11021) : the third equality also uses 
(11031) and the fact that D commutes with the diagonal matrix M 33 . Using the relation 
v_i = M_i it/i and the prescribed values of A±i, one may extract B±\. For example, if 
A_ x = 0, then 

Si(w) = C Ji .(a;)Ai(w) where C F (u) = HK e ™/ ys _ \ ■ ( 105 ) 

A rephrasing of the result (11051) is that the wave-function ty F (ui,t,y) suitable for computing 
the time-ordered propagator G F (u) = G u (w) has leading behavior 

G F (u) = [1 + C F (uj)y 3 }e~ iuJt (106) 
near the boundary of the R region. Then, in precise analogy to (131 

G F (u) = -g y/3 d log^ F (uj,t,y) 

' (107) 



-YiiXcot^ + 2 C ^ /yS - C ^ 



27 



where Y is defined as in ( I6T1) . Comparing with (1601) . one sees immediately that 

G F (u) = Re G R (u) + z coth — Im G R (u) . (108) 

tys 

Restoring factors of z H amounts to replacing ix by 1/T in (1 108j) . so this is the desired 
modification of the first equation in (1921) . 

Now let's return to the justification of the formulas (I99l - (ll03p . The main issue is contin- 
uing from ^±3 to \& ± 3 and from to ^I/.g. An argument was made in [53], and extended 
to the present context in [45], that the correct way to do this is to continue in the lower 
half of the complex V plane and the upper half of the complex U plane when horizon singu- 
larities are encountered. By so doing one imposes infalling conditions on positive frequency 
modes and outgoing conditions on negative frequency modes. Formally one may simply send 
U — > U + i5 and V — > V — i5 for small positive 5. In R this corresponds to sending y — > y + i5, 
and in L it corresponds to sending y — > y — iS. Consider the consequences for matching the 
oscillatory solution from 3 to 3, which takes place wholly within the R region: 

lim B 3 (y s -y- iSf^* = lim B~ 3 (y -y s + iS)™' 2 ** , (109) 
y->ys- y~*ys+ 

where the first limit is from below and the second is from above. From (11091) one obtains 
B3 = e -7rtJ / 2?/s I?3, whence the lower right entry in M33. The opposite behavior arises in 
M_ 3 _g because y — >• y — iS in the L region. In crossing the V = line, one finds (V — 
= e -W2(_y + j ( j)-<w/2. This completes a justification of ffTUD]) and (|I02jl . The 
other steps can be justified more straightforwardly, as follows: 

• f|99|) is a consequence of the wave equation for \l/ being the same in regions ±2 and ±4. 
The non-zero imaginary part of t doesn't affect the equation of motion. 

• (IIOH) is a consequence of applying the rule x ia = e lotl ° sx only when x is positive. The 
details are somewhat intricate, so I will trace the steps for the A-type solution from 
-5 to -5: 

A_ 5 (l - y)- <w / 4 e - <hrt =A- 5 {1- y )-W4 e -^Ret e -™/2 

/ ^2\^/ 4 (HO) 

where in the second equality I used the leading behavior of V near y — 1 in the L 
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region: 

{-Vf = (1 - y)e^ = (1 - y)e^^- . (Ill) 

The first and last expressions in fillip are in a form suitable to be raised to the imag- 
inary power —icu/A. 

(11031) is a consequence of the behavior of purely infalling and purely outgoing solutions 
of the wave equation for An infalling solution is one with B3 = 0, and its behavior 
the R boundary, (128jl . corresponds to V\ given by the first column of the matrix M13. 
An outgoing solution should have behavior at the R boundary corresponding to v% 
given by the second column of M43. It so happens that if ^>n(u,t,y) = e~ iut ij}ji(u},y) 
is the infalling wave-function constructed in section [3], then 

* A (u, t, y) = e-^f(y)^ 2 f(y/y s )-^ /2ys M^, vT (H2) 
is a purely outgoing solution, where 

f(y) = i±y e - 2 ^- 1 y. (113) 



Expanding (11121) for small y gives 

* A (u, t, y) = [l+ (C R (co)* - tX)y 3 ]e-^ , (114) 

as desired. 



C Speed limits on single quarks 

The string theory calculations presented in this paper pertain to external quarks: pointlike 
objects in the fundamental representation which have infinite mass associated with their 
near-field Coulombic color-electric field. Comparisons with QCD are thus better justified for 
c and b quarks, whose mass is well above the typical temperature of the QGP, than for light 
quarks. 

If the velocity of an external quark is too close to unity, there is a new reason to be 
wary of trailing string computations: the horizon on the worldsheet is very close to the 
boundary, where in a holographic representation of real QCD (assuming there is one) space 
presumably becomes highly curved, and the simplest calculations based on the Nambu-Goto 
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action in AdS'5-Schwarzschild may experience significant corrections. Let's try to estimate 
when problems of this sort might start to arise. It was suggested already in [24] to associate 
an energy scale 

_ L 2 /z H 1 _ T y/^VV (n5) 
27ra' y 2 y 

with a radius y in AdSs-Schwarzschild. This is justified by observing that a static string 
dangling from y into the horizon has mass 

^static = \ —T ( ~ ~ 1 ) ~ V for ?/ < 1. (116) 



2ira' \y 

Now suppose one identifies a scale /if a n at which AdS / CFT techniques (at least those based 
on supergravity and classical strings) start to fail. From (11151) one extracts a corresponding 
2/faii, and if the worldsheet horizon has ys < yMi there may be significant corrections to the 
trailing string results. Setting yg = yfaii and using (|23|) . one finds that the Lorentz factor of 
the heavy quark that the trailing string purports to describe is 

7M = AWvM • (117) 

The trouble with (jl 17j) is that gy M N, /x fail and T all incorporate considerable uncertainties. 

To get an idea of the range of plausible values for 7^, let's consider the two prescriptions 
f!74|) and (JTSj). I don't have a systematic way of determining /u,f a ii, but perhaps a reasonable 
range to consider is /if ai i = 0.8 — 1.6 GeV. Taking /i fai i = 1.2 GeV as a representative value, 
one finds 

"obvious:" 7M = 4.9( T |g v ) 



2 



"alternative:" 7fail = 29 (77^) 



(118) 



Because of the quadratic dependence on both /if ai i and T (not to mention the choice of 
comparison scheme) 7f a n remains substantially uncertain. Evidently, if the "obvious" pre- 
scription is used, there is some doubt cast on stringy predictions for charm quark when 
p c ^ 6.7GeV/c, corresponding to ~ 3.4 GeV. In the "alternative" scheme, there is less 
reason to worry about stringy corrections near the upper end of the trailing string. 

Another thing can go wrong if one wants to represent a finite mass quark as a string 
ending on a D7-brane, as in [4*8l [T]. Given m c for the charm quark, one can use (11161) to 
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obtain a corresponding position 

1 

Vc ~ 1 + 2m c /Ty/& M N (119) 

of the D7-brane. If ys < y c , then there is no horizon on the trailing string: the boundary 
of the worldsheet is spacelike. I regard this as a pathology which probably invalidates the 
trailing string picture for Lorentz factors 7 > 7 C = l/j/J?. Estimates of 7 C suffer from the 
same ambiguities as 7f a n, as discussed above. Using m c = 1.4 GeV and either (JTij) or (175]) . 
one arrives at the estimates 

"obvious:" 7 C = 13 

(120) 

"alternative:" 7 C = 53 . 

These values should again be regarded as incorporating considerable uncertainties: for ex- 
ample, it is puzzling that the horizon causes the mass of a quark to decrease from its zero- 
temperature value, but this decrease is what makes the values in (11201) higher than those 
derived from (11181) with /Xf a ;i — > m c . The main lesson to draw from fl 1 2 [) is that charm's 
mass is high enough to avoid threatening the existence of the worldsheet horizon for the 
momenta accessible at RHIC; but attempts to treat light quarks in terms of the trailing 
string construction are perilous indeed^ 



9 For example, in the "alternative" comparison scheme, an up quark whose mass in the medium is assumed 
to be m = 300 McV would be described in terms of a D7-brane with y u = 0.43, so ~f u = 5.5 is the maximum 
Lorentz factor. It is hardly fair to ignore corrections to relativistic dispersion relations in this context, but 
if one does so, the result is a total energy of 1.7 GeV. As an upper limit on allowed energies for partons 
described by trailing strings, this is fairly anemic. 



31 



References 



[1] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz, and L. G. Yaffe, "Energy loss of a 
heavy quark moving through N = 4 supersymmetric Yang-Mills plasma," 



hep-th/0605158 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 



[9 
[10 

[11 
[12 

[13 

[14 



S. S. Gubser, "Drag force in AdS/CFT," |hep-th/0605182 



J. Casalderrey-Solana and D. Teaney, "Heavy quark diffusion in strongly coupled N 



4 Yang Mills," hep-ph/0605199 



S.-J. Sin and I. Zahed, "Holography of radiation and jet quenching," Phys. Lett. B608 



(2005) 265-273, [hip^th/0407215 



H. Liu, K. Rajagopal, and U. A. Wiedemann, "Calculating the jet quenching 



parameter from AdS/CFT," hep-ph/0605178 



M. Chernicoff, J. A. Garcia, and A. Guijosa, "The energy of a moving quark-antiquark 



pair in an N = 4 SYM plasma," JEEP 09 (2006) 068, |hep-th/0607089 



H. Liu, K. Rajagopal, and U. A. Wiedemann, "Wilson loops in heavy ion collisions 



and their calculation in AdS/CFT," hep-ph/0612168 



P. C. Argyres, M. Edalati, and J. F. Vazquez-Poritz, "Spacelike strings and jet 



quenching from a Wilson loop," |hep-th/0 6 12157 



A. Buchel, "On jet quenching parameters in strongly coupled non-conformal gauge 



theories," Phys. Rev. D74 (2006) 046006, hep-th/0605178 



E. Nakano, S. Teraguchi, and W.-Y. Wen, "Drag force, jet quenching, and 



AdS/QCD," hep-ph/0608274 



P. Talavera, "Drag force in a string model dual to large-N QCD," hep-th/0610179 



Y.-h. Gao, W.-s. Xu, and D.-f. Zeng, "Jet quenching parameters of Sakai-Sugimoto 



model," hep-th/06 11217 



C. P. Herzog, "Energy loss of heavy quarks from asymptotically AdS geometries," 



JEEP 09 (2006) 032,|hep-th/0605191 



E. Caceres and A. Guijosa, "Drag force in charged N = 4 SYM plasma," JEEP 11 
(2006) 077, hep-th/0605235] 



32 



[15] F.-L. Lin and T. Matsuo, "Jet quenching parameter in medium with chemical 



potential from AdS/CFT," Phys. Lett. B641 (2006) 45-49, |hep^th70606136 



[16] S. D. Avramis and K. Sfetsos, "Supergravity and the jet quenching parameter in the 



presence of R-charge densities," hep-th/0606190 



[17] N. Armesto, J. D. Edelstein, and J. Mas, "Jet quenching at finite 't Hooft coupling 
and chemical potential from AdS/CFT," JEEP 09 (2006) 039, |hep-ph/0606245[ 

[18] E. Caceres and A. Guijosa, "On drag forces and jet quenching in strongly coupled 



plasmas," hep-th/0606134 



[19] J. F. Vazquez-Poritz, "Enhancing the jet quenching parameter from marginal 



deformations," hep-th/0605296 



[20] J. J. Friess, S. S. Gubser, and G. Michalogiorgakis, "Dissipation from a heavy quark 
moving through N = 4 super- Yang-Mills plasma," JEEP 09 (2006) 072, 



hep-th/0605292. 



[21] J. J. Friess, S. S. Gubser, G. Michalogiorgakis, and S. S. Pufu, "The stress tensor of a 



quark moving through N = 4 thermal plasma," hep-th/0607022 



[22] S.-J. Sin and I. Zahed, "Ampere's law and energy loss in AdS/CFT duality, 



hep-ph/0606049 



[23] M. Chernicoff and A. Guijosa, "Energy loss of gluons, baryons and k-quarks in an N 



4 SYM plasma," hep-th/0611155 



[24] J. M. Maldacena, "The large N limit of superconformal field theories and 



supergravity," Adv. Theor. Math. Phys. 2 (1998) 231-252, hep-th/9711200 



[25] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, "Gauge theory correlators from 



non-critical string theory," Phys. Lett. B428 (1998) 105-114, hep-th/9802109 



[26] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 (1998) 



253-291, hep-th/9802150 



[27] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri, and Y. Oz, "Large N field 



theories, string theory and gravity," Phys. Rept. 323 (2000) 183-386, hep-th/9905111 



33 



[28] E. D'Hoker and D. Z. Freedman, "Supersymmetric gauge theories and the AdS/CFT 



correspondence," hep-th/0201253 



[29] I. R. Klebanov, "TASI lectures: Introduction to the AdS/CFT correspondence," 



|hep-th/0009139 



[30] S.-J. Rey and J.-T. Yee, "Macroscopic strings as heavy quarks in large N gauge theory 



and anti-de Sitter supergravity," Eur. Phys. J. C22 (2001) 379-394, [hep-th/9803001 
[31] J. M. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80 (1998) 



4859-4862, hep-th/9803002 



[32] STAR Collaboration, B. I. Abelev et. at, "Transverse momentum and centrality 
dependence of high-p(T) non-photonic electron suppression in Au + Au collisions at 
v /s^7 = 200 GeV," |nucl-ex/0607012| 

[33] PHENIX Collaboration, A. Adare et. ai, "Energy Loss and Flow of Heavy Quarks in 



Au+Au Collisions at ^/s~^ = 200 GeV," |nucl-ex/0611018 



[34] S. S. Gubser, "Comparing the drag force on heavy quarks in N = 4 super- Yang-Mills 



theory and QCD," hep-th/06 11272 



[35] H. van Hees, V. Greco, and R. Rapp, "Heavy-quark probes of the quark-gluon plasma 
at RHIC," Phys. Rev. C73 (2006) 034913, jnHTl-th/05080551 

[36] B. Svetitsky, "Diffusion of charmed quark in the quark-gluon plasma," Phys. Rev. 
D37 (1988) 2484-2491. 

[37] H. van Hees and R. Rapp, "Thermalization of heavy quarks in the quark-gluon 
plasma," Phys. Rev. C71 (2005) 034907, lmTTl-th/04120151 



[38] G. D. Moore and D. Teaney, "How much do heavy quarks thermalize in a heavy ion 



collision?," Phys. Rev. C71 (2005) 064904, [56^11/0412346 



R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne, and D. Schiff, "Radiative energy 
loss of high energy quarks and gluons in a finite- volume quark-gluon plasma," Nucl. 



Phys. B483 (1997) 291-320, hep-ph/9607355 



[40] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne, and D. Schiff, "Radiative energy 
loss and p(T)-broadening of high energy partons in nuclei," Nucl. Phys. B484 (1997) 



265-282, hep-ph/9608322l 



31 



[41] B. G. Zakharov, "Radiative energy loss of high energy quarks in finite-size nuclear 



matter and quark-gluon plasma," JETP Lett. 65 (1997) 615-620, |hep :: ph/9704255 



[42] K. Peeters, J. Sonnenschein, and M. Zamaklar, "Holographic melting and related 
properties of mesons in a quark gluon plasma," Phys. Rev. D74 (2006) 106008, 



hep-th/0606195 



[43] H. Liu, K. Rajagopal, and U. A. Wiedemann, "An AdS/CFT calculation of screening 



in a hot wind," hep-ph/0607062 



[44] D. T. Son and A. O. Starinets, "Minkowski-space correlators in AdS/CFT 



correspondence: Recipe and applications," JEEP 09 (2002) 042, hep-th/0205051 



[45] J. Casalderrey-Solana and D. Teaney, "Transverse Momentum Broadening of a Fast 



Quark in a N=4 Yang Mills Plasma," hep-th/0701123 



[46] C. G. Callan and A. Guijosa, "Undulating strings and gauge theory waves," Nucl. 
Phys. B565 (2000) 157-175, [hip-th/9906153| 

[47] A. Dainese, C. Loizides, and G. Paic, "Leading-particle suppression in high energy 



nucleus nucleus collisions," Eur. Phys. J. C38 (2005) 461-474, hep-ph/0406201 



[48] A. Karch and E. Katz, "Adding flavor to AdS/CFT," JEEP 06 (2002) 043, 



hep-th/0205236 



[49] S. J. Brodsky and P. Hoyer, "A Bound on the energy loss of partons in nuclei," Phys. 



Lett. B298 (1993) 165-170, hep-ph/9210262 



[50] K. Schweda et. al, "D anti-D correlations as a sensitive probe for thermalization in 



high-energy nuclear collisions," |nucl-ex/ 0610 043 , 

[51] N. Armesto, M. Cacciari, A. Dainese, C. A. Salgado, and U. A. Wiedemann, "How 
sensitive are high-p(T) electron spectra at RHIC to heavy quark energy loss?," Phys. 



Lett. B637 (2006) 362-366, hep-ph/0511257 



[52] PHENIX Collaboration, A. Adare et. al, "Measurement of high-p(T) single electrons 



from heavy-flavor decays in p + p collisions at y/s = 200 GeV," hep-ex/0609010 



[53] M. Cacciari, P. Nason, and R. Vogt, "QCD predictions for charm and bottom 



production at RHIC," Phys. Rev. Lett. 95 (2005) 122001, hep-ph/0502203| 



35 



[54] C. P. Herzog and D. T. Son, "Schwinger-Keldysh propagators from AdS/CFT 



correspondence," JEEP 03 (2003) 046, |hep-th702T2072 



[55] G. T. Horowitz and D. Marolf, "A new approach to string cosmology," JEEP 07 



(1998) 014, |hep-th/9805207 



[56] J. M. Maldacena, "Eternal black holes in Anti-de-Sitter," JEEP 04 (2003) 021, 



hep-th/0106112 



[57] V. Balasubramanian, P. Kraus, A. E. Lawrence, and S. P. Trivedi, "Holographic probes 



of anti-de Sitter space-times," Phys. Rev. D59 (1999) 104021, [5^^/9808017 



36 



